An Eccentric Hot Jupiter Orbiting the Subgiant HD 185269 1 

John Asher Johnson 2 , Geoffrey W. Marcy 2 , Debra A. Fischer 3 , Gregory W. Henry 4 , Jason 
T. Wright 2 , Howard Isaacson 3 , Chris McCarthy 3 



O 
O 

(N 

> 
O 

> 
in 
m 
O 
oo 
O 

o 



9 L,: 
6 

CO 

C3 



j ohnj ohnOastron . berkeley . edu 



ABSTRACT 



We report the detection of a Jupiter-mass planet in a 6.838 d orbit around 
the 1.28 M subgiant HD 185269. The eccentricity of HD 185269b (e = 0.30) is 
unusually large compared to other planets within 0.1 AU of their stars. Pho- 
tometric observations demonstrate that the star is constant to ±0.0001 mag on 
the radial velocity period, strengthening our interpretation of a planetary com- 
panion. This planet was detected as part of our radial velocity survey of evolved 
stars located on the subgiant branch of the H R diagram — also known as the 
Hertzsprung Gap. These stars, which have masses between 1.2 and 2.5 M , play 
an important role in the investigation of the frequency of extrasolar planets as a 
function of stellar mass. 

Subject headings: techniques: radial velocities — planetary systems: formation- 
stars: individual (HD 189269) 



X 



1. Introduction 

The assemblage of ~ 185 known extrasolar planets^ has revealed important relation- 
ships between the physical characteri stics of stars and the likelihood that they harbor planets 



flUdrv et all 12003 



Butler et al.l 120061). One example is the strong correlation between plane t 



occurrence and stellar metallicity ( Gonzalezlll997l ; ISantos et all2004j : lFischer fc Valenti 



2005J). 



1 Bascd on observations obtained at the Lick Observatory, which is operated by the University of California 
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1 For the updated catalog of extrasolar planet and their parameters see http://exoplanets.org 
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The planet-metallicity relationship can be understood in the context of the core accretion 
model of planet formation, in which Jovian planets begin as large rocky cores that then ac- 
crete gas_^^m_Jh^irsuiTounding protoplanetary disks once they exceed a critical core mass 



(e.g. iPollack et al.lll996l ). The growth of these embryonic cores is enhanced by increasing 
the surface densit y of solid particles in the disk, whic h is related to the metallicity of the 
star/disk system (jAlibert et al.ll2005l ; llda fc Linll2005al ). 



The growth rate of rocky cores in protoplanetary disks is also related to the total disk 
mass. Assuming that the disk mass increases with the mass of the central star, the planet 
occurrence rate should therefore also correlate with stellar mass. lLaughlin et al.l (120041 ) 
showed that the lower surface densities of M dwarf protoplanetary disks impede the growth 
of Jupiter-mass planets. T he relationship betw een stellar mass and planet occurrence was 
studied in further detail by llda fc Lin! (l2005bl ) for a larger range of stellar masses. Based 
on their Monte Carlo simulations, they predict a positive correlation between the number of 
detectable planets and stellar mass up to about 1 M Q . However, there is some debate about 
whether surface density of solid material in protoplanetary disks is proportional to the mass 
of the central star . By assuming that the initial conditions of the disk are independent from 
the mass of star, iKornet et al.l (120061 ) find that the occurrence rate of planets is inversely 
related to stellar mass. 

The prediction that low-mass stars should have a lower frequency of Jovian planets is in 
accordan ce with observation s; only one M dwarf, GL 876, is known to harbor Jupiter -mass 
planets (IMarcy et al.l 120011 ). From their survey of 90 M dwarfs, lEndl et al.l (120061 ) esti- 
mate that l e ss tha n 1.27% of such low-mass stars harbor Jovian-mass planets within 1 AU. 
Laws et al.l (120031 ) studied the planet occurrence rate in light of the current observational 
data using the larger set of FGK stars surveyed as part of the California fe Carneg ie Planet 
Search (CCPS, the full, updated target list can be found in IWright et al.l (120041 ))) . They 
found evidence that the planet rate decreases for lower stellar masses and peaks near 1.0 M G . 
However, the planet rate is poorly constrained for stellar masses greater than about 1.2 M Q 
due to the small number of intermediate-mass stars (1.3 < M star < 3 M ) in the CCPS 
sample. This dearth of massive stars is due to an observation al bias since main-sequenc e 
stars with spectral types earlier than F8 tend to be fast rotators (Ido Nascimento et al.112 003') 



have fewer spectral lines, and display a large amount of chromospheric activity (jSaar et al. 



19981 ). These features result in a decrease in the radial velocity precision attainable from the 
spectra of main-sequence stars more massive than ~ 1.2 M jGalland et al.l 120051 ). Thus, 



early-type dwarfs are not typically monitored as part of most radial velocity surveys. 



One method to circumvent these difficulties is to observe intermediate-mass stars after 
they evolve into the region of the H-R diagram between the main-sequence and red giant 
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branch — also known as the Hertzsprung Gap (HG). After stars have expended their core 
hydrogen fuel sources their radii expand, their photospheres cool, and convection sets in 
below their photospheres. Convective motion and stellar rotation drive magnetic fields that 



couple with an expanding stellar wind and act a s a rotational brake (IGray fc Nagarlll985 



Schrijver fc Polall993l ; Ido Nascimento et al.ll2000l ). The cooler atmospheres and slower rota- 
tional velocities of evolved stars lead to an increased number of narrow absorption lines in 
their spectra, making HG stars better suited for precise radial velocity measurements than 
their main-sequence progenitors. 

We are conducting a radial velocity survey of 159 HG stars to search for planets orbiting 
intermediate-mass stars. We describe the selection criteria of our sample of stars in § |2j We 
present here the first planet detection from our sample of HG stars: an eccentric hot jupiter 
orbiting the 1.28 M subgiant, HD 185269. The properties of the host star are presented in 
§ |3l and we describe our observations and orbital solution in § HI We present the results of 
our photometric monitoring in § [5] and conclude with a discussion in § [6l 



2. Sample 



We have selected 159 HG stars based on the crite ria 0.5 < M y < 3.5, 0.55 < B-V < 1.0, 
and V < 7.6, as listed in the Hipparcos Catalog (lESAl 119971 ). Additionally, we selected 
only stars lyin g more than 1 mag above the mean Hipparcos main-sequence, as defined by 
Wrightl (120041 ) . We chose the red cutoff to avoid red gi ants, which are known to exhibit 
excess velocity jitter (IFrink et al.ll2001uHekker et al.ll2006l ). The lower My restriction avoids 
Cepheid variables, and the upper limit excludes stars with masses less than 1.2 M Q . We 
avoided stars in the clump region of the H-R diagram (B — V > 0.8 and My < 2.0) in 
order to avoid the mass ambiguity stemming from the closely-spaced, overlapping isochrones 
between low-mass horizontal branch stars and high-mass stars on their first ascent to the 
giant branch. 
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Fig. 1. — Our sample of Hertzsprung Gap stars (filled circles) with respect to the mean 
Hipparcos main-sequence (thick line) and the theoretical mass tracks of iGirardi et al.l (120021 ) 
for [Fe/H] = 0.0. The mass tracks for Solar-mass stars on the horizontal branch with 
[Fe/H] = {—0.7, —0.4,0.0} (from left to right) are shown as dashed lines in the upper right 
portion of the figure. The overlapping mass tracks in this region of the H-R diagram leads 
to ambiguious mass estimates. The filled star shows the position of HD 185269. The hashed 
regions show the approximate ranges of the samples of other planet search programs. 
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The H-R diagram of our full sample is illustrated in Figured], along with the mass tracks 
of iGirardi et al.l (120021 ) . The mass range of our sample is 1.2 < M* < 2.5M , with a mean 
of 1.5 M . We are observing 115 of the brightest and northernmost of these stars at Lick 
Observatory and the remaining 44 at Keck Observatory. Figure [TJ also shows the approxi- 
mate search domains of other programs containing evolved, intermediate-mass stars. These 



surve ys include se arches for planets around clump gian ts (ISato et al.l 120031 ; ISetiawan et al. 



20031 ). red giants (IFrink et al.l 120021 ; lHatzes et al.l 120051 ) and the subgiants included in the 



CCPS. As evidenced from Figure [TJ HG stars occupy a unique and unexplored region in the 
H-R diagram. 



3. Properties of HD 185269 



HD 185269 (=HIP 96507) is a G0IV subgiant with V = 6.67, B — V = 0.60 6, a parallax- 
based distance of 47.6 pc, and an absolute magnitude My = 3.29 (lESAlll997l ). Its position 



in the H-R diagram in relation to other stars in our sample is shown in Figure [TJ HD 185269 
l ies 1.1 mag a bove the mean Hipparcos main-sequence of stars in the solar neighborhood 
( iWrightl 120041 ). confirming its subgiant classification. This star is chromospherically-quiet 



with S = 0.14 and R' HK = —5.14, from which we estimat e a 23 d rotation period. We used 



the LTE spectral synthesis described by lValenti &: Fischer! (120051 ) to calculate T e ff = 5980 K, 
[Fe/H] = +0.11, logg = 3.94 and V rot smi = 6.1 km s -1 . We interpol ated the star' s color , 
absolute magnitude and metallicity onto t he stellar model gri ds of IGirardi et al.l (120021 ) 
using the Bayesian methodology detailed by iPont fc Eyerl (120041 ) . Our interpolation yields 



an estimated stellar mass M* = 1.28 ± 0.1 M Q , radi us R* = 1.88 ± 0.07 in agree ment 



with M* = 1.31 M and R* = 1.74 R Q estimated bv lAUende Prieto fc Lambertl (I1999I ). All 
of the stellar properties are summarized in Table [TJ 



4. Observations and Orbit 

We began observing HD 185269 in 2004 May at Lick Observatory using the 0.6 m Coude 
Auxiliary Telescope (CAT) and the 3 m Shane telescope. Both telescopes feed the Hamilton 
echelle spectrograph, thus no offset correction is needed for the velocities. We measured 
radial veloci ti es fro m the high-resolution (X/5X = 50,000) spectral observations using the 



Butler et all (119961 ) iodine cell method. Traditionally, this method requires an additional 
reference observation, or template, made without the iodine cell. These template observa- 
tions require higher signal and resolution than normal radial velocity observations, which 
leads to increased exposure times. Given our large target list and the small aperture of the 
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Fig. 2. — Phased radial velocity measurements of HD 185269. The dashed line shows the 
best-fit orbital solution. The bottom panel shows the phased residuals of the fit, which have 
RMS = 10.1 m s" 1 . 

CAT, obtaining an observed template for each star would represent a prohibitive cost in 
observing time. We instead analyze the initial observations of a ll our targets u s ing sy nthetic 



(or "morphed") templates, following the method described by iJohnson et al.l (120061 ). Stars 
showing excess RMS scatter are reanalyzed using a high-quality observed template, obtained 
with the 3 m Shane telescope, to achieve improved long-term velocity precision. 

The first five observations of HD 185269 spanned 1 yr and had an RMS scatter of 
53 m s -1 , prompting us to initiate follow-up observations with an increased sampling rate 
to search for a short-period signal. Our radial velocity measurements are listed in Table [2] 
along with the Julian dates and internal measurement uncertainties. A periodogram analysis 
of the velocities reveals a 6.838 periodicity with a false alarm probability < 0.001%. To 
search for a full orbital solution, we augmented our measurement uncertainties with a jitter 
estimate of 5.0 m s _1 , based on the star's chromospheric activity index, absolute magnitude 
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and color (ISaar et al.l Il998l ; IWrightl 120051 ). The best-fit Keplerian orbital solution yields 
an orbital period P = 6.838 d, velocity semi-amplitude K = 91 m s _1 and eccentricity 
e = 0.30 ± 0.04. Using our stellar mass estimate of 1.28 M Q , we derive M P sin i = 0.94 M Jup 
and a = 0.077 AU. Figure [2] shows our velocities and orbital solution phased with the 6.838 d 
period. The full orbital solution and parameter uncertainties are listed in Table [31 

We estimate the parameter uncertainties using a Monte Carlo method. For each of 
100 trials the best-fit Keplerian is subtracted from the measured velocities. The residuals 
are then scrambled and added back to the original measurements, and a new set of orbital 
parameters is obtained. The standard deviations of the parameters derived from all trials 
are adopted as the la uncertainties listed in Table [31 



Photometric Observations 



Queloz et al.l (120011 ) and iPaulson et al.l (120041 ) have shown that active regions such as 



spots and plages on the photospheres of solar-type stars can cause low-amplitude, radial ve- 
locity variations (jitter) by distorting the stellar line profiles as the spots are carried across 
the stellar disk by rotation. If a large active region exists for several stellar rotations (not 
an unusual circumstance for stars younger and more active than the Sun), then periodic 
rotational modulation of the spectral line profiles can mimic the presence of a planetary 
companion. Therefore, precision photometric measurements can be an important comple- 
ment to Doppler observations. For radial velocity variations cause d by surface magn etic 
activity, the star will exhibit low-level photometric variability (e.g., iHenry et al.l 119951 ) on 
the radial velocity period. If the radial velocity variability is the result of true reflex motion 
caused by a planetary companion in orbit around the star , the star in genera l will not show 



photometric variability on the radial velocity period (e.g.. IHenry et al.ll2000al ). Photometric 
observations of planetary-candidate stars can also detect transits of planetary companions 
with inclinations near 90° and so allow the determination of a planet's true mass, rad ius, 



density, and composition (e.g.. IHenry et al.ll2000bl ; ISato et al.ll2005l ; iBouchy et al.l 120051 ) 
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Planetary Orbital Phase 

Fig. 3. — Stromgren (b + y)/2 photometric observations of HD 185269 acquired with the 
Til 0.8 m APT at Fairborn Observatory. The observations have been phased to the orbital 
period of the planet and an estimated time of mid transit. {Top) There is no evidence 
for any periodicity in the observations between 1 and 100 days. In particular, the star is 
constant on the radial velocity period to a limit of 0.0001 mag or so, supporting the planetary 
interpretation of the radial velocity variations. (Bottom) The photometric observations 
around the predicted time of transit replotted with an expanded scale on the abscissa. The 
predicted transit depth of only 0.003 mag is shown for an assumed planetary radius of 1 
Rj up ; such transits are ruled out by the observations. 



- 9- 



We observed HD 185269 with the Til 0.8 m automatic photometric telescope (APT) 
at Fairborn Observatory during May and June of 2006, obtaining a total of 207 brightness 
measurements. The Til APT is equipped with a two-channel precision photometer employ- 
ing two EMI 9124QB bi- alkali photomultiplier tubes to make simultaneous measurements in 
the Stromgren b and y passbands. The APT measures the difference in brightness between 
a program star and a nearby constant comparison star or stars with a typical precision of 
0.0015 mag for bright stars (V < 8.0). For HD 185269, we used the two comparison stars 
HD 184151 (V = 6.87, B - V = 0.46, F5 V) and HD 184381 (V = 6.70, B - V = 0.45, 
F5 V). Differential magnitudes between the two comparison stars showed them both to be 
constant to 0.0012 mag or better on a night-to-night timescale. We created Stromgren b and 
y differential magnitudes of HD 185269 with respect to the average of the two comparison 
stars to improve our photometric precision. The differential magnitudes were reduced with 
nightly extinction coefficients and transformed to the Stromgren system with yearly mean 
transformation coefficients. To improve precision still further, we combined the separate b 
and y differential magnitudes into a single (b + y)/2 pass band. Additional information on 



Henrv 


(1999) and 


Eaton et al. 


(2003) 



The 207 combined (b + y)/2 differential magnitudes of HD 185269 are plotted in the top 
panel of Figure [31 The observations are phased with the planetary orbital period and the 
time of mid transit given in Table [3l The standard deviation of the observations from the 
mean brightness level is 0.0011 mag, suggesting that HD 185269 as well as its comparison 
stars are all highly constant. Period analysis does not reveal any periodicity between 0.03 
and 100 days. A least-squares sine fit of the observations phased to the radial velocity 
period gives a semi- amplitude of 0.00015 ± 0.00012 mag. This very low limit to possible 
photometric variability supports planetary-reflex motion as the cause of the radial velocity 
variations. 

In the bottom panel of Figure [3], the observations near phase 0.0 are replotted with an 
expanded scale on the abscissa. The solid curve in each of the two panels approximates the 
predicted transit light curve assuming a planetary orbital inclination of 90° (central transits). 
The out-of-transit light level corresponds to the mean brightness of the observations. The 
predicted transit duration is calculated from the orbital elements, while the predicted transit 
depth of 0.003 mag is derived from the stellar radius of 1.88 R from Tableland an assumed 
planetary radius of 1.0 Rj up . Thus, any transits of the planet across the subgia nt star are 



expec ted to be very shallow, as was the case with the transits of HD 149026 (ISato et al. 



20051 ). The horizontal bar below the predicted transit window represents the approximate 
uncertainty in the time of mid transit, computed from the orbital elements. The vertical 
error bar to the right of the transit window corresponds to the ± 0.0011 mag measurement 
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uncertainties for a single observation. The geometric probability of transit s in this system 



i s ~1 2%, computed from the orbital elements in Tableland equation 1 of ISeagroves et al. 



( 120031 ). Although the uncertainty in the time of mid transit is slightly larger than the 
predicted duration of any transits, the observations nonetheless rule out the existence of 
transits except perhaps for short events occurring around phase 0.97. In the absence of 
transits, the orbital inclination must be less than ~ 83°. 



6. Summary and Discussion 



We are monitoring the radial velocities of a sample of 159 intermediate-mass Hertzsprung 
Gap (HG) stars in order to study the relationship between stellar mass and planet occurrence 
rate. We present the detection of a 0.94 Mj up planet in an eccentric, 6.838 d orbit around 
the 1.28 M subgiant HD 185269. 

Compared to other planets with orbital separations a < 0.1 AU (also known as "hot 
jupiters"), HD 185269b has an unusually large eccentricity. Figure [H shows th e distribution 
of eccentricities for the 33 planets having a < 0.1 AU (IButler et al.l 120061 ). With e = 
0.30, HD 185269 b stands out from the distribution as one of only three hot j upiters with 
eccen tricities in excess of 0.2. The o ther two planets are HP 162020 b (e = 0.277; lUdry et al. 
2002h and HD 118203 b (e = 0.309: Ida Silva et alJbood ). 



The origin of eccentricities amongst single exoplanets is not well understood, but the 
gravitational influence of additional companions can drive the eccentricities of planets in mul- 
ticomponent systems. An as-yet undetected sec ond planet is likely the cause of the large ec- 



centricity of HD 118203b — Ida Silva et al.l (120061 ) report a linear trend with a 49.7 m s 1 yr 



slope. However, we do not see evidence of a linear trend in our observations of HD 185269, 
nor significant periodicities in the Keplerian fit residuals. We are continuing to monitor this 
star at Keck and Lick Observatories to search for additional low-mass companions. 

It is thought that the nearly circular orbits of most hot jupiters is a result of tidal 
circularization. This leaves open the possibility that the tidal circularization timescale at 
the orbital distance of HD 185269b is longer than the age of the star. Indeed, theoretical 
predictions of the circ ularization timescale have a strong dependence with semi-major axis, 
typically t circ oc a 13//3 (ITerquem et al.lll998l ). Eight of the 10 planets with a > 0.055 AU have 
eccentricities greater than 0.1, which may indicate that the timescale for tidal circularization 
at distances beyond ~ 0.055 AU is longer than the age of most FGK s tars. The detection 
of additional hot jupiters by programs such as the N2K Consortium (IFischer et al.l 120051 ) 
will help shed light on the relationship between the eccentricity and orbital separation of 
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Fig. 4. — The eccentricity-period distribution of hot jupiters listed in iButler et al.l (120061 ). 
HD 185269 is shown as a filled pentagram. The inset shows the histogram of eccentricities 
for the same sample of stars. The arrow denotes the position of HD 185269. 



short-period planets. 

Other sea rches for planets around intermediate -mass stars have so far focused on the red 



giant branch (IFrink et al. 



flSato et alJl200a 



200 



2: 



Setiawan et al 



20o3). These programs have to date discovered a total of 



iatzes et al.l 120051 ) and clump regions of the H-R diagram 



6 sub-stellar objects orbiting giant stars, proving that planets do form and can be detected 
around evolved stars more massive than ~ 1.5 M Q . However, stars in the clump and red giant 
branches follow closely-spaced and often overlapping evolutionary tracks, making precise 
stellar mass estimations from isochrone interpolation difficult (see Figured]). By contrast, 
HG stars have mass tracks that are nearly parallel and widely-spaced in My, enabling precise 
mass deter minations. HG sta r s also exhibit lower velocity jitter and have smaller radii than 
red giants (IFrink et al.l l200ll ; iHekker et al.l 120061 ) , which enables the detection of a wider 
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range of planet masses and orbital separations. We will present the stellar characteristics 
and velocity behavior of our sample in a forthcoming paper. 
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including Shannon Patel, Julia Kregenow, Karin Sandstrom, Katie Peek and Bernie Walp. 
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tion of the Lick Observatory staff. We appreciate funding from NASA grant NNG05GK92G 
(to GWM) for supporting this research. DAF is a Cottrell Science Scholar of Research 
Corporation and acknowledges support from NASA Grant NNG05G164G that made this 
work possible. GWH acknowledges support from NASA grant NCC5-511 and NSF grant 
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Table 1. Stellar Parameters 



V 


6.67 


My 


3.29 (0.08) 


B-V 


0.606 


Spectral Type 


GOIV 


Distance (pc) 


47 (1.0) 


[Fe/H] 


+0.11 (0.05) 


Teff (K) 


5980 (50) 


V rot sini (km s _1 ) 


6.1 (0.5) 


\ogg 


3.94 (0.07) 


■U. (M ) 


1.28 (0.1) 


-R* (-R©) 


1.88 (0.1) 




0.14 




-5.14 


Prot (d) 


23 


Age (Gyr) 


4.2 


Jitter (m s _1 ) 


5.0 
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Table 2. Radial Velocities for HD 185269 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



13155.980 


-6.5 


5.8 


13211.839 


30.7 


6.2 


13522.875 


-107.0 


5.6 


13576.840 


-34.9 


7.2 


13619.734 


-100.0 


5.7 


13629.661 


51.7 


6.2 


13668.621 


-21.5 


6.7 


13718.641 


49.6 


6.3 


13720.587 


-34.3 


6.9 


13866.932 


-23.5 


6.4 


13867.878 


11.8 


4.3 


13868.806 


38.9 


5.9 


13868.996 


44.3 


5.6 


13869.982 


32.4 


6.8 


13879.841 


-63.7 


5.4 


13880.912 


16.0 


7.1 


13884.918 


-50.7 


5.6 


13885.924 


-148.0 


5.5 


13891.795 


-66.5 


6.9 


13891.963 


-85.7 


7.0 


13894.908 


32.8 


4.9 


13895.895 


44.6 


5.0 


13896.895 


42.2 


4.5 


13897.908 


-13.8 


5.1 


13898.905 


-78.0 


12.0 


13900.929 


-31.9 


5.2 


13901.947 


26.7 


5.7 


13902.915 


40.5 


4.6 


13903.936 


29.1 


5.6 


13904.888 


-25.6 


6.0 



Table 2 — Continued 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



Table 3. Orbital Parameters for HD 185269 b 



P(d) 


6.838 (0.001) 


V (JD) 


2453154.1 (0.18) 


T transit (JD) 


2453153.1 (0.16) 


e 


0.30 (0.04) 


Ki (m s- 1 ) 


91 (4.5) 


uj (deg) 


173 (6.8) 


M P sin i (M Jup ) 


0.94 


a (AU) 


0.077 


Fit RMS (m s- 1 ) 


10.1 


A / 

V ^-v 


1.14 


N obs 


30 



a Time of periastron passage. 



